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Abstract 

Recent numerical work on the Zabusky-Kruskal experiment has revealed, amongst other things, the existence of 
hidden solitons in the wave profile. Here, using Osborne's nonlinear Fourier analysis, which is based on the periodic, 
inverse scattering transform, the hidden soliton hypothesis is corroborated, and the exact number of solitons, their 
amplitudes and their reference level is computed. Other "less nonlinear" oscillation modes, which are not solitons, 
are also found to have nontrivial energy contributions over certain ranges of the dispersion parameter In addition, 
the reference level is found to be a non-monotone function of the dispersion parameter. Finally, in the case of large 
dispersion, we show that the one-term nonlinear Fourier series yields a very accurate approximate solution in terms 
of Jacobian elliptic functions. 
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1. Introduction 

Four decades after the discovery of solitons by Zabusky and Kmskal (ZK) ll35l through a computational experi- 
ment, the study of the evolution of harmonic initial data under the Korteweg-de Vries (KdV) equation on a periodic 
interval is far from complete [4|. Beyond the discovery |35 1 of the ability of the KdV equation's (localized) traveling- 
wave solutions (termed 'solitons') to retain their "identity" (shape, speed) after collisions, modern numerical simula- 
tions by Salupere et al. of this paradigm equation of solitonics reveal such exotic features as "hidden" (or "virtual") 
solitons li9J, emergence of soliton ensembles li32J and long-time periodic patterns of the trajectories ||29l[33l . These 
phenomena have been shown to be generic of nonUnear waves, as they occur under other governing equations as well 

This raises the simple, yet quite fundamental, question: How many solitons emerge from a harmonic input? A 
successful approach to answering this question is based upon discrete spectral analysis ifTOl [30l |3T1 . The essence of 
this method is to characterize the solitary waves based on the information inherent in the pseudospectral numerical 
approximation of the underlying partial differential equation (PDE) |26|. In general, solitary wave identification is 
a difficult problem ifTSl [36l . especially when multiple wave interactions occur and long time scales are considered. 
Fortunately, in the case of the ZK experiment, the KdV equation has an advantage over other nonlinear wave equations 
in that it is integrable, i.e., it can be solved exactly (in theory) using the inverse scattering transform UJ on both the 
infinite line and on a periodic interval. 

In this respect, Osborne's nonlinear Fourier analysis 1 16] provides a natural framework for applying the periodic, 
inverse scattering transform (PIST) for the KdV equation to real-world problems. In particular, it overcomes the 
difficulty of the PIST being only a theoretical tool by providing a practical numerical implementation of it lfT7l[T9ll24l . 
Osborne and Bergamasco ETIl employed this approach to successfully reproduce the numerical results of Zabusky and 
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Kruskal 1351 . They were able to confirm the number of solitons observed emerging from a harmonic input and the 
recurrence time of the initial condition. In present work, we employ the same approach to corroborate the numerical 
results of Salupere et al. regarding hidden solitons. Specifically, the aim here is to determine the hidden modes' 
amplitudes and classify them in the hierarchy of solutions to the periodic KdV equation (i.e., as either solitons, 
cnoidal waves or harmonic waves). Conceptually, one may consider this approach as a generalization of some of the 
ideas in aT0ll3T1 . 

This paper is organized as follows. In Sec. [2] the physical form of the Korteweg-de Vries equation and its 
transformation to the form considered in |35 | is presented. In Sec. [3] the interpretation of the PIST as nonlinear 
Fourier analysis is discussed. In Sec. [4j the PIST spectrum of the harmonic initial condition for the KdV equation is 
shown for various values of the dispersion parameter, and the hidden soliton hypothesis is discussed. Sec. |5] in the 
large-dispersion case, illustrates in more detail how a nonlinear Fourier series is constructed using the PIST. Finally, 
before concluding in Sec.|7] in Sec.|6]we elaborate on the notion of a soliton reference level for the periodic problem 
and its dependence on the dispersion parameter 



2. Position of tlie problem 

The KdV equation is arguably the most famous "soliton-bearing" PDEs, governing phenomena as seemingly 
disparate as the motion of lattices, the collective behavior of plasmas and the shape of hydrodynamic waves (see, e.g., 
Ifn |7] [22] [35j and the references therein). The classical context 1 14| in which it arises is the propagation of "long" 
waves over "shallow" water. (For purposes of the present work, we do not need to make these terms any sharper.) 
Letting rjix, t) be the surface elevation, the KdV equation in the moving frame is 

rjt + Ci)!], + aT]T],- + = 0, (x, t) e [0, L] x (0, oo), (1) 

where L{> 0) is the length of the domain, and the subscripts denote partial differentiation with respect to an indepen- 
dent variable. In addition, the speed of linear waves cq, the nonlinearity coefficient a and the dispersion coefficient fi 
are constant physical parameters. For surface water waves, for example, they can be expressed in terms of the channel 
depth h and the acceleration due to gravity g as follows ll22l : 

ry 3co cqIi^ 

CO = ^Jgh, 

As in ll35l . we are interested in the periodic Cauchy problem (i.e., the initial-value problem subject to periodic 
boundary conditions). Thus, we have that 

ri{x,0) ^ t]q{x), x€[0,L], 
T](x + L,t) ^ Tj(x,t), (x,f) e [0,L] X [0,oo). 

A particularly illustrative initial condition, which is the one we focus on in the present study, is the harmonic one: 

2n 

Tio(x) - acos(a)x), tu - — n, n e Z, (4) 

L 

where a is an arbitrary (real) constant. 

Now, while Eq. ([T]) is in the appropriate form to apply Osborne's nonlinear Fourier analysis, it is not in the 
form originally considered by Zabusky and Kruskal 1351 and the subsequent works of Salupere et al. To make the 
comparison possible, we introduce the following new variables: 

1 T 

T](x, t) - -u{^, t), x-(^ + Cot) mod L, t - t, (3 - 6 , (5) 

a 

and for definiteness we take L - 2 cm and n - \. Upon substituting the latter transformations into Eq. ([TJ, we obtain 

Uj + uu£ + 5^u^^^^Q, (^,t) e [0,2] X (0,oo), (6) 
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and the initial-boundary conditions in Eq. Q become 



m(^, 0) = aa cosing), ^ e [0, 2] , 

+ 2, t) = m(^, t), , t) e [0, 2] x [0, oo). ^'^^ 



Clearly, we must choose a - l/a to normalize the initial condition as in | 

Also, we must rewrite the parameters in Eq. ([!} in terms of the known quantities 6 and g cm/s^. To this end, 

we use the identity fi - 6^ and Eq. ([2]) to deduce 

^ (8) 



Finally, we note that the dispersion parameter 6 used here is identical to the one in ll35l . and it is related to the 
dispersion parameter d\ of Salupere et al. via di - -2 logig(7r(5). 



3. Interpretation of the PIST as a "nonlinear Fourier transform" 

Next, we turn to the relationship between the PIST and the (ordinary) Fourier transform, and the interpretation 
of the former as a nonlinear generalization of the latter. To this end, first we note that the solution strategy by the 
scattering transform can be split into two distinct steps: the direct problem and the inverse problem. The former 
consists of solving the eigenvalue problem 

•H(A = £i/', 9{ -.^-—r + Vix), xe[0,L], (9) 

where V(x) := -Arjoix) is the "potential," A - a/(6yS) is a nonlinearity-to-dispersion ratio, and £ e M is a spectral 
eigenvalue. Equation (|9| has been studied extensively: in quantum mechanics it is the celebrated (time-independent) 
Schrodinger equation fTI|, and, in the theory of ODEs, it is known as Hill's equation f5l. For periodic "potentials," 
i.e., when V(x + L) - V{x) Vx e [0, L] as we have assumed, it is well-known that the spectrum of the operator 'H 
is divided into two distinct sets depending upon the boundary conditions imposed on the eigenfunctions ij/ ||T1|8| [T91 . 
Thus, it is common to classify the spectral eigenvalues (also known as the "scattering data") as belonging to either the 
main spectrum, which we write as the set {£j)^^j^', or the auxiliary spectrum, which we write as the set {ju'^j^j, where 
is the number of degrees of freedom (i.e., oscillations modes IIT6I or band gaps ID). 
On the other hand, the inverse problem consists of constructing the nonlinear Fourier series from the spectrum 
U [jl'j] using either Abelian hyperelliptic functions |18, 23| or the Riemann ©-function |3, 20|. In former case, 
which is the so-called yu-representation of the PIST, the exact solution of Eq. Q, subject to the initial and boundary 
conditions given in Eq. ([3]), takes the form 

^ ( N 2N+1 \ 

riix, = - 1 2 g ^ij{x, t)-Yj&j\. (10) 

It is important to note that all nonlinear waves and their nonlinear interactions are accounted for in this linear 
superposition. Unfortunately, the computation of the nonlinear oscillation modes (i.e., the hyperelliptic functions 
{pj{x,t)]^^y) is highly nontrivial; however, numerical approaches have been developed \Q3\ and successfully used in 
practice fa |25]| . 

Several special cases of Eq. ( [TO| i off'er insight into why the latter is analogous to the ordinary Fourier series 
and aid with the interpretation of the results in the following sections. In the small-amplitude limit, i.e., when 
n\aXxj\nj{x,t)\ « 1, we have fJ.j(x,t) ~ cos(kjX - ojjt + (pj), where kj is the wavenumber, is the frequency 
and 4>j the phase of the mode. Therefore, if we suppose that all the oscillation modes fall in the small-amplitude 
limit, then Eq. ( fTO] ) reduces to the ordinary Fourier series! This relationship is more than just an analogy, a rig- 
orous derivation of the (ordinary) Fourier transform from the scattering transform, in the small-amplitude limit, is 
given in |17|. Next, if there are no interactions, e.g., the spectrum consists of a single wave (i.e., = 1), we have 
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fiiix, f) - cn^{k\x - (jL>\t + (pi\mi), which is a Jacobian elliptic function with modulus mi. In fact, it is the well-known 
cnoidal wave solution of the (periodic) KdV equation |1 1. 

For the hyperelliptic representation of the nonlinear Fourier series, given by Eq. ([T0|, the wavenumbers are com- 
mensurable with those of the ordinary Fourier series, i.e., kj - Inj/L (I < j < N) IfTSll . However, this is not the only 
way to classify the nonlinear oscillations. One can use the modulus nij, termed the "soliton index," of each of the 
hyperelliptic functions, which can be computed from the main spectrum as 

- - ^'^-'-^'^ l<j<N. (11) 



£>2j+l - &2j-\ 

Then, each nonlinear oscillation mode can be placed into one of two distinct categories based on its soliton index: 

1. Mj > 0.99 solitons, in particular, cn^(f |ot = 1) = sech^(^); 

2. ntj «: 1.0 => radiation, in particular, cv?{^\m = 0) = cos^(^). 

For Mj not in either distinguished limit above, the qualitative structure of the nonlinear oscillations is not immediately 
obvious. Boyd [T| argues that for moderate moduli the polycnoidal wave solutions of the KdV equation are actually 
well-approximated by the first few terms in their Fourier series, showing they are, in fact, not much different from the 
linear waves of the mj <sc 1 .0 limit. 

Furthermore, it can be shown Il22ll25l that the amplitudes of the hyperelliptic functions are given by 



-(Ei-s;f - £>2j), for solitons; 

\ (12) 
— i^2j+i - &2j), Otherwise (radiation). 



where iiref = ^f+i is the soliton reference level with /* being the largest j for which mj > 0.99. Then, clearly, the 
number of solitons in the spectrum is A^soi = ./*■ 

Finally, we note that the theory of Eq. (|9]l is quite mature and exact solutions can be obtained for a number of 
specific forms of the potential V{x) |1J|5|. Unfortunately, for the ZK experiment (recall Eq. Q), we have V(x) - 
-/(acos(wx) for which there is no known closed form solution. All is not lost, however, the essence of Osborne's 
nonlinear Fourier analysis is that a numerical solution of Eq. (|9| can be obtained and the details of the PIST carried 
out in this way. To this end, we use a modified version of Osborne's automatic algorithm 1 19], as described in |6 1. The 
result is an exact (to any desired numerical precision) representation of the nonlinear Fourier spectrum of any initial 
condition of the KdV equation. 



4. How many solitons in a cosine wave? 

4.1. The Zabusky-Kruskal experiment 

Recall that, in Sec.|2]upon switching to the new set of variables, we distilled all physical parameters of the KdV 
equation into the dispersion parameter 6. Therefore, it suffices to vary S to establish all possible ways a harmonic 
initial condition can evolve under the KdV equation. In this subsection, we analyze the classical ZK experiment, i.e., 
6 - 0.022, using the PIST. To this end, in Fig. [T| the ordinary and nonlinear Fourier spectra of the harmonic initial 
condition are presented. Since the ordinary Fourier transform, implemented here using the Fast Fourier Transform 
(EFT), does not resolve the temporal evolution of a solution of the KdV equation, the EFT of the solution, which was 
computed by direct numerical integration of the PDF, is presented at f = 3.6/n. This temporal value was chosen for 
ease of comparison with [35] and because all observable solitons are visible in the wave profile at this instant of time. 

The EFT spectrum suggests that there are over thirty significant normal modes of the problem, with the majority 
of the enegry concentrated in three distinct wave number bands centered at k x 2.5 cm"'. A: ^ 25 cm"' and k x 47.5 
cm"'. Of course, without further study of the EFT spectrum and its evolution in time, nothing can be said about the 
number of solitons present in the solution profile. 

In contrast, the PIST shows exactly eight soliton modes plus other nonlinear waves and radiation. As was shown 
in on . the PIST spectrum captures precisely the eight solitons observed in |35|. However, what has never been 
discussed before are the other (four of them, in fact, as Table [T] below shows) nontrivial modes in the spectrum. 
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FFT Spectrum at f = 3. 6(3 = 3.61k 
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Figure 1 : (Color online.) Comparison of the ordinary Fourier spectium (left panel) at a specific time and the nonlinear Fourier spectrum (right 
panel) of the original ZK experiment {S = 0.022). The vertical dashed line in the right panel denotes the wave number of the last true soliton mode. 

By 'nontrivial' we mean their amplitudes are not so small as to have negligible energy contribution. The discrete 
spectral analysis approach of Salupere et al. ||30l l3T| found there exist at least four hidden "solitons" beyond the eight 
observable ones. Indeed, in Fig.[T] an additional four modes beyond the "true" soUton ones are easily distinguished 
(i.e., have large enough amplitude). Thus, the hidden "soliton" hypothesis is true. What is more, the PIST naturally 
classifies these hidden modes within the hierarchy of solutions to the periodic KdV equation; it happens they are not 
proper solitons (i.e., they have mj < 0.99). 

Finally, Fig. [T] illustrates a very important point: namely, that the FFT is only useful when analyzing a solution. 
That is to say, the underlying PDE has to be solved up to some instant of time, and then the spectrum computed from 
this data. Other tools, that can provide both time and frequency resolution, such as the short-time (or windowed) 
Fourier transform exist fV2\, however, the PDE must still be solved in advance. On the other hand, the PIST fully 
characterizes the initial condition and its evolution. This is due to the fact that the PIST is a formal technique for 
integrating exactly the KdV and other such integrable nonlinear wave equations. 

4.2. The effect of dispersion on soliton generation 

It has been shown ifTOl [30l |3TI that the number of solitons detected in the solution of the KdV equation varies 
with the dispersion parameter Indeed, it is well-known that, as the zero-dispersion limit (6 — » 0) is approached, the 
number of solitons generated by an initial condition grows quickly [, 1 1 J . In Fig.|2j we present the PIST spectra of the 
solution for a number of representative values of 6 besides the ZK value; these were chosen so that a comparison with 
ifTOl 1301 13 II is possible. The quantitative results from theses figures, and all other ones presented in the his paper, are 
also summarized in Table[T] which includes the relevant data from lfT0l[30l[3Tl and the heuristic estimate A^soi ~ 0.2/6 
duetoZabuskvlMl. 

For 6 - 0.0178999 (o dj - 2.5), we observe eleven solitons and nine non-soliton nonlinear waves for a total of 
twenty waves generated by the harmonic initial condition. Though this value of 6 was too small for the computational 
efforts at the time of the studies in ifTOllSTl . the PIST has no issues whatsoever with 5 — > 0. The computational time 
for the nonlinear Fourier analysis grows only linearly with the number of spectral eigenvalues, and it is independent 
of 6. 

For 6 = 0.0252843 (o di - 2.2), we observe six solitons and eight non-soliton nonlinear waves. For 6 - 
0.0635112 (o di = 1.4), the PIST finds two solitons and six other modes. And, for 6 = 0.142184 (o di = 0.7), 
there are zero solitons out of a total of five waves in the spectrum. The pattern is clear: as the dispersion parameter 
increases the eflfects of nonlinearly weaken and fewer solitons are produced from the initial harmonic wave. The 
number of other non-soliton modes in the spectrum also decreases with 6, from nine at 5 = 0.0178999 to five at 
6 = 0.142184. 

Now, returning to Table [T] and comparing the number of solitons Nso\ given by the PIST and the number of visible 
solitons observed by Salupere et al. (i.e., total minus hidden), we see that all but one of the observed "solitons" in the 
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S = 0.0178999 <5 = 0.0252843 
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Figure 2: (Color online.) Nonlinear Fourier spectra of the harmonic initial condition for four representative values of 6. As before, a vertical dashed 
line delineates the last soliton mode from the remainder of the spectrum. 

wave profile are indeed solitons, and the last is a highly-nonlinear cnoidal wave. This shows that across a wide range 
of 6 values, as was the case for the ZK experiment discussed in the previous subsection, the hidden "solitons" are a 
manifestation of the other nonlinearly interacting non-soliton solutions of the periodic KdV equation. 
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Table 1 : Comparison of the discrete spectral and nonlinear Fourier analyses of hidden "solitons" and extended statistics of the nonlinear Fourier 
spectrum for all values of 5 considered here. 



Finally, we note that one cannot expect to detect (in, e.g., a numerical experiment) all nonlinear oscillation modes 
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the PIST finds. That is, the discrete spectral analysis ifTOl [3TI is inherently limited in its ability to resolve low am- 
plitude, or very weakly nonlinear waves in the spectrum. Therefore, to shed some light on precisely which modes 
can be detected and why, the last two column of Table |T| show the number NAj>m-^ of modes with amplitudes greater 
than 10"^ and the number N„y^io-\ of modes with elliptic moduli greater than 10 respectively. Though these cut- 
offs are largely arbitrary, they are qualitatively reasonable. Indeed, the total number of "solitons" found by Salupere 
et al. is either of these numbers, for the values of 6 considered here. Thus, the "small" and "weak" modes in the 
spectrum (typically, these are radiation modes) are difficult to resolve through discrete spectral analysis. However, we 
have shown, beyond any doubt, that the discrete spectral analysis properly distinguishes the total number of highly 
nonlinear oscillation modes, both soliton and otherwise, and that there indeed are hidden modes in the ZK experiment. 



5. A note on the large dispersion case 

Though, from a mathematical point of view, the zero-dispersion (6 — > 0) limit is by far the most interesting distin- 
guished limit of the KdV equation IfTTII . the large-dispersion one (6 = (9(1)) illustrates very well some of the concepts 
behind the PIST and Osborne's nonlinear Fourier analysis. To this end, Fig.|3]shows the nonlinear Fourier spectrum 
of the harmonic initial conditions for two large values of 6. Clearly, as 5 — > 1, there is only one large-amplitude mode, 
with the remaining two modes' amplitudes approaching zero as 6 increases (recall Table[T]l. Moreover, the "degree of 
nonlinearity" (i.e., the elliptic modulus nij) of the leading mode decreases as 6 increases, while its amplitude remains 
the same. This is due to the fact that the terms 6^u^^^ and uu^ in Eq. (|6]l now balance asymptotically, and there is no 
nonlinear steepening of the wave profile that leads to the generation of sohtons ll35l . 



(5 = 0.317999 6 = 1. 
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Figure 3: (Color online.) Nonlinear' Fourier spectra of the harmonic initial condition for large values of the dispersion parameter The vertical 
dashed line representing the end of the soliton spectrum is now aXk = because solitons do not emerge in the solution for these values of S. 



Thus, to a good approximation, N = I when 6 = 0(1). In this case, the solution of the periodic KdV equation 
takes the following form IfTSl : 

u{^,t) X 4Ai cn^iki^ + LJiT\mi) + A-\6 - 263), (13) 

where ki - yj&T, - &\ is the mode's wavenumber, ux - 2&ki is its frequency, and we have set & :- Yj^i £>j for 
convenience. Of course, {£;}J^p A\ - 0.499987 and mi - 0.0653120 are obtained from the PIST. This approximate 
solution is compared to the numerical solution of the KdV equation, for (5 = 1, in Fig.|4] The agreement between 
the two is quite good. More importantly, Eq. ( [TJI illustrates a fundamental difference between the infinite-line and 
periodic versions of the 1ST. Not only does the periodic problem have a much richer solution space (a continuum of 
solutions expressible as Jacobian elliptic functions), but it also requires three spectral eigenvalues to determine one 
mode. This means that the infinite-line 1ST approach in IfTOllSTI . while helpful, cannot capture the full picture of the 
periodic problem. 
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Figure 4: (Color online.) Comparison between the numerical (blue) and approximate PIST solution (dashed green), given by Eq. \13) , of the 
large-dispersion KdV equation (r5 = 1). 

6. A reference level for the periodic problem 

Recall that in Eq. ( [T2] i, we defined the quantity E^gf as the soliton reference level. This idea was first formulated 
in on . where Osborne and Bergamasco showed that a shift in the reference (or zero) level of a multi-soliton solution 
is mathematically equivalent to a shift in the energy level of the last soliton's band gap edge with respect to £ = 0. 
This gives the natural definition = £2/+! ~ 0, where (as in Sec.[3]l / is the largest j for which mj> 0.99. Hence, 
if we know E^^i, then we can determine the solitons' reference propagation level and reference wavenumber for the 
periodic problem as follow: 

"ref = -aE,^flA, = 2nj* jL. (14) 

Table|2]gives Mret as a function of the dispersion parameter for all values of 6 discussed in the preceding sections. 
Naturally, in the large-dispersion case no solitons form, and the harmonic initial condition's shape is largely preserved 
as it propagates. In this limit, the propagation level is most naturally understood as the mean level of the solution. 
Indeed, it is very close to zero — the mean of a trigonometric function over its period. This leads to an important 
observation: in the presence of large radiation modes in the spectrum, the reference level is not simply the minimum 
of the wave profile in the solution. That is to say, the true Wief cannot be observed in the numerical solution of the 
problem, as the oscillatory radiation modes will continually shift the wave profile minimum; if they could somehow 
be "filtered" one would see the "pure" solitons propagating on the level given by Mret in Table |2] 
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Table 2: The soliton reference level H,ef , computed from the PIST spectrum using Eq. for various values of the dispersion parameter 5. 

What is more interesting, however, is that Mref increases with 6, until at a certain value it reverses itself and 
goes back through zero becoming negative. As 5 ^ 0, it appears that the PIST predicts that Mref — > -1, which is 
the "absolute minimum" level observed numerically f TOl l30l |3TI . This shows that the PIST and direct numerical 
simulation approaches are consistent, as they should, but that the PIST off'ers a deeper insight into how Mref varies with 
6, which is quite difficult to infer from the numerical simulation. 

Finally, we note that the infinite-line inverse scattering transform (1ST) was used in IfTOl 1311 to corroborate the 
discrete spectral analysis (numerical) approach. But, because the periodic problem for the KdV equation is not 
mathematically equivalent to the infinite-line problem with an initial condition of compact support, the choice of Mref 
in lfT0l [3T1 had to be made before the number of solitons could be calculated. This lead to an artificial dependence 
of A^soi on Mref. Though good results were obtained using an appropriate estimate for Mref, it is clear that the PIST 
approach, which is the natural one for the periodic problem, gives an unambiguous way to calculate Mref. 
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7. Conclusion 



The present work shows that the evolution ("spectrum") of harmonic initial data under the Korteweg-de Vries 
(KdV) equation can be fully and automatically characterized by the periodic, inverse scattering transform (PIST). In 
particular, within the framework of Osborne and Bergamasco |21 1, the hidden soliton hypothesis of Salupere et al. is 
corroborated, and the exact number of solitons, their amplitudes and their reference level is computed. This offers new 
insight into the phenomenon because such precise results were not possible through the direct numerical integration 
of the KdV equation and the discrete spectral analysis approach [101 123 111] [30l |3T1 [32l |33] . In particular, the 
apparent linear variation of the soliton amplitudes with the wavenumber, first discussed in ll35l . is true for all soliton 
modes detected by the PIST, over a range of 6 values. Meanwhile, the remainder of the spectrum consists of other 
"less nonlinear" oscillation modes, e.g., nonlinearly-interacting cnoidal waves. Though the amplitudes of these do not 
follow the linear trend, they are not negligible, and they are, in fact, what is observed numerically as "hidden solitons." 

It is important to note that the present approach (via the PIST and Osborne's nonlinear Fourier analysis) to the 
problem of soliton formation is complementary to the discrete spectral analysis approach of Salupere et al. While 
the PIST allows for a very precise decomposition of the initial condition into the "basis" elements of the nonlinear 
PDE at hand (in the present case, the KdV equation), it requires that the PDE be fully integrable — a rather stringent 
stipulation. The discrete spectral analysis approach to soliton formation does not require any a priori mathematical 
structure of the PDE. 

Finally, two of the findings of the present work require further investigation. First, the non-monotone dependence 
of Mief on 6 is quite unexpected. A detailed study is necessary to identify the mechanism for this. Second, the number 
A^soi of solitons predicted by the PIST is always less than the number of "solitons" found in the numerical solution by 
Salupere et al. across the entire range of values of 6 considered here. This means that most hidden "solitons" are not 
solitons, per se. Clearly, the richness of the solution space of the periodic KdV equation allows for modes that fall 
"in-between" solitons and radiation. Whether this distinction can be made using the discrete spectral analysis is a very 
interesting question that would be very relevant when studying non-integrable wave equations, for which a nonlinear 
Fourier transform cannot be constructed, and therefore the distinction between soliton and non-soliton modes cannot 
be made analytically. 
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